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ABSTRACT 

We have combined deep photometry in the B,V and I bands from CTIO/MOSAIC of the Sculptor dwarf 
spheroidal galaxy, going down to the oldest Main Sequence Turn-Offs, with spectroscopic metallicity 
distributions of Red Giant Branch stars. This allows us to obtain the most detailed and complete Star 
Formation History to date, as well as an accurate timescale for chemical enrichment. 
The Star Formation History shows that Sculptor is dominated by old (>10 Gyr), metal-poor stars, but that 
younger, more metal-rich populations are also present. Using Star Formation Histories determined at dif- 
ferent radii from the centre we show that Sculptor formed stars with an increasing central concentration 
with time. The old, metal-poor populations are present at all radii, while more metal-rich, younger stars 
are more centrally concentrated. We find that within an elliptical radius of 1 degree, or 1.5 kpc from the 
centre, a total mass in stars of 7.8xl0 6 M Q was formed, between 14 and 7 Gyr ago, with a peak at 13-14 
Gyr ago. 

We use the detailed Star Formation History to determine age estimates for individual Red Giant Branch 
stars with high resolution spectroscopic abundances. Thus, for the first time, we can directly determine 
detailed timescales for the evolution of individual chemical elements. We find that the trends in alpha- 
elements match what is expected from an extended, relatively uninterrupted period of star formation con- 
tinuing for 6-7 Gyr. The knee in the alpha-element distribution occurs at an age of 10.9± lGyr, suggesting 
that SNe la enrichment began «2±lGyr after the start of star formation in Sculptor. 

Key words. Galaxies: dwarf - Galaxies: evolution - Galaxies: stellar content - Galaxies: Local Group - 
Stars: C-M diagrams 
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1. Introduction 

The Sculptor dwarf spheroidal galaxy (dSph) is a 
relatively faint (My«— 11.2), well studied system 
in the Local Group. A distance of 86+5 kpc has 
been determined from RR Lyrae star measurements, 
in good agreement with other distance determina- 
tions, such as t he tip of the RGB an d horizon- 
tal branch level (Pietrzv fiski et alJ f2008. and refer- 
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ences therein). Sculptor is located at high galac- 
tic latitude (b=- 83°) with relatively l ow reddening, 
EfB-V)=0.018 (ISchleeel et alJ[l998h . 
Past studies of the Star Formation History (SFH) 
of Sculptor have been made using colour-magnitude 
diagrams (CMDs) of varying depth, and spatial cov- 
erage. Using CMDs down to the Main Sequence 
Tur n-Off (MSTO) in a field just outside the core ra- 
dius iDa Costal dl98l determined an age range for 
Sculptor of 13+2 Gyr. A small field of view (~2'), 
deep HST image of a field well outside the cen- 
tre of Sculptor confirmed the ancient age (15+2 
Gyr) of the bulk of the stars in Sculptor using 
synthesis CMD analysis, but also showed a small 
tail of star formation reaching down to more re - 
cent times dMonkiewicz et al.lll999HDolphinl2d02h . 
An effort was made to combine detailed spectro- 
scopic abundances of 5 stars with photometric ages, 
which also suggested that Sculptor is predomi- 
nantly _oJd J _witha small tail of intermediate age 
stars dTolstov et alj|2003h . 

In a previous qualitative study of wide-field CMDs, 
the Horizontal Branch (HB) morph ology was found 
to change significantly with ra dius (Maiews ki et alj 
119991: lHurlev-KelleretalJll999h. which was late r 
linked to a metallicity gradient (Tols toy etalj2 004). 
From deep wide-field CMDs reaching the oldest 
MSTOs, covering a large fraction of Sculptor, this 
metallicity grad ient was also found to be linked to 
an age gradient dde Boer et alj|20l"it) . 
In Sculptor, wide-field medium resolution Can 
triplet spectr oscopy of a large number of RGB stars 
is available (Battagli jgoTfft iBattaglia etaT]|2008bt 
Starkenbur g et al.l2010l) . giving a well defined spec- 
troscopic Metallicity Distribution Function (MDF) 
for stars with ages >1.5 Gyr old. In the central 25' 
diameter region of Sculptor high res olution (HR) 
spect roscopy (Hill et al., in prep, see iTolstov et alj 
2009) of stars on the upper RGB provides detailed 
abundances of a-elements (O, Mg, Ca, Si, Ti) as 
well as r- and s-process elements (Y, La, Ba, Eu, 
Nd). Furthermore, within the central 15' diameter 
region stars have been observed using medium res- 
olution spectroscopy going dow n to fainter mag- 
nitudes dKirbv et alj|2009l I2010h . giving [Fe/H] as 
well as a-element abundances. 
The presence of these gradients suggests that 
Sculptor has experienced an initial, spatially ex- 
tended episode of metal-poor star formation at all 
radii, and stars of higher metallicity were subse- 
quently formed more towards the centre. This pic- 
ture explains the properties of the CMD, such as the 
different spatial distributions traced by the HB mor- 
phology and Red Giant Branch (RGB) metallicity 




Fig. 1: Coverage of the photometric and spectroscopic ob- 
servations across the Sculptor dwarf spheroidal galaxy. 
The solid (black) square denotes the full coverage of 
the CTIO 4m/MOSAIC fields. The (blue) open circles 
show the stars observe d in the VLT/FLAMES low resolu- 
tion Can triplet survey (Battaglia 2007; Starkenburg et al. 
120101) . The (red) solid dots mark the RGB stars that 
also have high reso lution abundance m easurements (Hill 
et al., in prep, see ITolstov et al] 2009) and the (green) 
da shed ellipse is the tidal radiu s of Sculptor, as given 
bv llrwin & Hatzid imitriov] il995h . 



gradient. 

Sculptor is a galaxy that can be used as a bench- 
mark for a simple galaxy with an extended episode 
of star formation in the early Universe. However, 
in order to quantify this picture and obtain accurate 
timescales for different stellar populations at differ- 
ent radii the detailed SFH needs to be determined 
over a large area of the galaxy. 
The Sculptor dSph has been modelled sev- 
eral times in simulations u s ing d i fferent tech- 
nique s (e.g.. [Salvadori et alj 120081: iRevaz eta! 
2009; Lanfran chi & Matteuccil 720101; iKirbv et al 
12011 : IReva z & Jabl onka 1201 lh . The cosmological 



semi-analitycal model of Salv adori et al.1 d2008h fol- 
lows simultaneously the evolution of the Milky 
Way and its dwarf galaxy satellites, and repro- 
duces the observed MD F and total mass con- 
tent of the Sculptor dSph. lLanfranchi & Matteuccil 
(2010) use chemical evolution modelling to re- 
produce the observed MDF of a variety of dwarf 
spheroidal galaxies, including Sculptor (s ee also 
lLanfranchi & Matteuccil 12004) . IKirbv eTal] (1201 lh 
use a chemical evolution model to match the ob- 
served MDF and alpha-element distributions from a 
large sample of spectroscopic ob servations in Local 
Group dwarf spheroidal galaxies. lRevaz & Jablonkal 
(1201 lh use a chemo-dynamical Smoothed-particle 
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Hydrodynamics (SPH) code to model the proper- 
ties of Local Group dwarf spheroidal galaxies. Their 
model for the Sculptor dSph correctly matches 
the observed MDF and shows an extended SFH. 
Furthermore, the narrow [Mg/Fe] distribution in 
their Sculptor-like model matches well with obser- 
vations from HR spectroscopy. 
In this paper the SFH of a large area of 
the Sculptor dSph will be d etermined using 
CMP synthesis methods (e.g., Tbsi et alJ 1 1 99 lb 



Tolsto y & Sahll996t|Gallart e t al. 1996a; Dolphin 
19971 : lAparicio etld] 1 19971) to inte rpret the deep 
wide-f ield photometry presented in Ide Boer et alJ 
d201 ll) . The available spectroscopic information 
from Can triplet and HR spectroscopy will be di- 
rectly used together with the photometry to provide 
additional constraints on the SFH. The spatial cov- 
erage of the photometric data (covering «80% of the 
tidal radius of Sculptor) allows us to determine the 
SFH at different radii and quantify the observed ra- 
dial age and metallicity gradients. 
Using the detailed SFH we also determine the prob- 
ability distribution function for age for stars on 
the upper RGB, giving age estimates for individ- 
ual stars. By linking these ages to the observed 
spectroscopic abundances we directly obtain, for the 
first time, the timescale of enrichment from different 
types of Supernovae (SNe). 

The paper is structured as follows: in section [2] we 
present our photometric and spectroscopic observa- 
tions. In section [3] we describe our method of ob- 
taining the SFH using the synthetic CMD method, 
adapted to include the MDF information. The de- 
tailed SFH analysis of the Sculptor dSph is given 
in section [4] The details of the age determination 
of individual stars are given in Section 14741 Finally, 
the conclusions that can be drawn from the SFH are 
discussed in section[5] 

2. Data 

Deep optical images of the Sculptor dSph in the B, 
V and I bands were obtained using the CTIO 4-m 
MOSAIC II camera. The reduction and accurate 
calibration of this d ataset is described in detail in 
a preceding paper dde Boer et alJ l20l"Tb . The total 
coverage of the nine photometric fields observed 
is shown in Figure Q] as the solid black square. 
The spatial coverage of the B band photometry is 
complete for radii r e //<0.5 degrees, while the V and 
I bands are complete for r e //<l degree. 
By stacking together several images for each 
pointing the deepest photometry possible was 
obtained. Short exposures were also obtained, to 



be able to accurately photometer the bright stars 
that are saturated in the deep images. In order to 
ensure accurate photometric calibration, images 
were obtained with the 0.9m CTIO telescope 
under photometric conditions. Observation s were 
also made of Landolt standard fields dLandoltl 
I2007L 119921) . Photometry was carried out on these 
images using DoPHOT dSchechter et alJ Il993h . 
The different fields were placed on the same 
photometric scale and combined in order to create a 
single, care fully calibrated pho tometric catalog, as 
described in de Boe r et al.l d2011l) . 




-1.7 -1.5 
[Fe/H] (dex) 

Fig. 2: Spectroscopic [Fe/H] and [a/Fe] measurements 
with their errors obtained from HR observations (Hill et 
al., in prep) in the central 25' diameter of the Sculptor 
dSph. The solid (black) line shows the range of possible 
values of [a/Fe] assumed for each bin of [Fe/H] in the 
SFH determination. 



In addition, Can triplet spectroscopy is available 
for »630 individual RGB stars that are likely mem- 
bers of Sculptor, from medium re solution (R~6500) 



VLT/FLAMES obser v ations dBatt aglia 



iBattaglia etal] l2008bt IStarkenburg et al 
These observations measure [Fe/H] for 



2007; 
l2010h . 



measure [J^e/HJ tor a large 
sample of stars for which we also have photometry, 
out to a radius of 1.7 degrees from the centre of the 
Sculptor dSph (see Figure[T). 
Furthermore, HR spectroscopy (VLT/FLAMES) 
is available for 89 individual RGB stars within 
the central part of the Sculptor dSp h, for r e /;<0.2 
degre es (Hill et al., in prep, see iTolstov et alJ 
2009). These observations provide [Fe/H] as well 
as [a/Fe] measurements. For [a/Fe] we assume 
[a/Fe] =( [Mg/Fe] +[Ca/Fe] +[Ti/Fe])/3. The HR 
spectroscopy includes a range in metallicity from 
-2.5<[Fe/H]<-1.0 dex, with alpha-element abun- 
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dances showing a clear correlation with [Fe/H] for 
the range -0.3<[a/Fe]<0.5 dex (See Figure[2]l. 



3. Method 

CMDs contain the signatures of numerous evolu- 
tionary parameters, such as age, chemical abun- 
dance, initial mass function, etc. The typical method 
of determining a SFH relies on comparing observed 
and synthetic CMDs of the individual stars that 
can be resolved. There have been many schemes 
proposed to quantify the SFH extracted from the 
CMP (e g.jTosii 



atity t 

et ajj|199lt iTolstov & Sahalll996h 
Gallart et all 1 1 996at iDolphinl 1 1 9971: 1 Aparicio et atl 



1997tlDolphirJ2002l) . 

We have created our own routine, Talos, which uses 
an approach that compares observed CMDs with a 
grid of synthetic CMDs through the use of Hess di- 
agrams (density plots of stars in the CMP) to deter- 
mine the SFH (similar to iDolphinll 1 997l 120021) . The 
synthetic CMDs include observational effects in a 
statistical manner, which provides the most realistic 
way to include them. 

We assume that a SFH can be built up from a lin- 
ear combination of simple stellar populations. The 
advantage of this approach is that instead of synthe- 
sising a large number of artificial CMDs, each with 
their own complex Star Formation Rate (SFR(t)), 
the simple populations only need to be generated 
once. Then, the combination of simple stellar pop- 
ulations has to be found, which best represents the 
observed CMD. 

Using the standard technique a number of different 
CMDs (e.g., V,B-V and I,V-I) can be independently 
used to obtain SFHs. A CMD is inherently two di- 
mensional (colour and magnitude), whereas when 
photometric information is available in more than 
two filters, we can use 3D information to more pre- 
cisely constrain the SFH. 

To fully utilize all the available photometric infor- 
mation to constrain the SFH, Talos fits the measure- 
ments in all the available passbands at once. In this 
way all the advantages of the different CMDs (such 
as the more precise photometry in the I,V-I CMD 
and the larger colour range of the V,B-I CMD) are 
incorporated into a single, more accurate SFH. 
In addition to the photometric data we can also add 
spectroscopic observations of a large number of in- 
dividual RGB stars. To take into account this extra 
information Talos also fits the observed MDF at the 
same time as the photometry, which allows us to put 
well motivated constraints on the metallicity range 
of the different stellar populations. 



The SFH determination technique consists of the 
following steps: 

1. Construct synthetic CMD and MDF models 

2. Find the combination of models that best match 
the data. 

3. Determine the resolution and uncertainties of 
the final SFH. 

These steps will be explained in more detail in the 
following sections. 

3.1. Constructing synthetic CMD models 

The synthetic stellar evolution library adopted 
in Talos is the D arthmouth Stellar Evolution 
Database dDotter et al.l2008l) . The database contains 
isochrones distributed with a specific grid of age, 
metallicity and a-element abundance, among which 
are also those with [a/Fe]<0. These are of particular 
importance for our analysis, since Sculptor contains 
a substantial number of stars (33% of the total) with 
low ff-element abundances (See Figure|2|. 
In order to allow greater flexibility in choosing the 
population domain and parameter resolution to use 
when modelling the SFH, we interpolate between 
the provided isochrones to produce a finer grid. A 
routine is provided by the Dartmouth group, to in- 
terpolate isochrones in [Fe/H]. Interpolation in age 
and [ff/Fe] is done linearly, given that the grid of 
age points is sufficiently fine that there will only be 
small changes between isochrones of different age. 
The isochrone library does not model the HB, 
Asymptotic Giant Branch (AGB) or Blue Straggler 
Star (BSS) phases, which means we cannot use 
these evolutionary features to constrain the SFH. 
The AGB is especially important, since it merges 
with the RGB in the CMD. The largest uncertain- 
ties in our final model may come from misiden- 
tifying the AGB using RGB models or the BSS 
as a young population. However, the Dartmouth 
isochrones have been shown to give a good si- 
multaneous fit to all other evolutiona ry features 
within a CMD, justifying their use (e.g jGlatt et alj 
l200&3lbh . 

Synthetic CMDs are generated by drawing stars 
from the isochrones according to an initial 
mass function (IMF) and the mass range within 
each isochrone. T he IMF used is the Kroupa 
IMF dKroupal2001l) . An initial SFR for a stellar pop- 
ulation (within a default mass range of 0. 1-120 M ) 
is assumed (resulting in several times more syn- 
thetic stars than in the observed CMD), which en- 
sures that enough stars are generated in each part of 
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the synthetic CMD, consistent with the time spent 
in each evolutionary phase. These synthetic stars are 
then placed at the distance of the Sculptor dSph, as- 
suming a reddening of E(B-V)=0.018. 
This gives us ideal CMDs for each stellar population 
that we consider in the SFH analysis. To create a 
synthetic CMD that can be compared directly to the 
observed photometry we need to include observa- 
tional effects in order to simulate our observational 
limitations. 

3.1.1. Adding observational effects 

The crucial aspect in comparing observed and syn- 
thetic CMDs is to determine the observational bi- 
ases, such as photometric errors, incompleteness, 
etc. 

This is done by carrying out a large number of sim- 
ulations, in which a number of artificial stars (with 
known brightness) are placed on the observed im- 
ages. These images are then re-reduced in exactly 
the same way as the original images, after which 
the artificial stars are recovered from the photom- 
etry. In this way we obtain a look-up table, which 
can be used to accurately model the effects of ob- 
servational condition s (e.g JStetson & Harrisll 1988k 
iGallart et alj Il99 6b). The lookup-table is used in 
Talos to assign an individual artificial star (with 
similar colours and magnitudes) to each star in an 
ideal synthetic CMD and considering the manner in 
which this star is recovered to be representative of 
the effect of the observational biases. 




1.2 1.4 1.6 



Fig. 3: The recovered I,V-I CMD of the artificial stars put 
into the observed images of Sculptor. The colours indicate 
the completeness fraction of artificial stars in each bin, 
with a scale on the right hand side of the plot. 



Using individual artificial stars has the disadvantage 
that a huge number of stars need to be generated in 
order to correctly sample the effects in each CMD 



bin. This is because using the same artificial star 
in assigning offsets to several synthetic stars cre- 
ates artificial clustering in the CMD. However, this 
approach is the only way to properly take into ac- 
count the colour-dependence of the completeness 
level and the asymmetry of the offset s applied to 
model stars at faint magnitudes (e.g. IGallart et al.l 
Il996bh . 
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Fig. 4: The magnitude of the 50% completeness level with 
increasing elliptical radius (r>//) from the centre of the 
Sculptor dSph, in B (solid blue line), V (green dashed 
line) and I (red dotted line) filters. The vertical (black) line 
indicates the elliptical radius up to which observations in 
the B filter are available. 



The set of artificial stars was generated with pa- 
rameters encompassing the range 5<Age<15 Gyr, 
-2.5<[Fe/H]<-0.80 dex, -0.2<[a/Fe]<0.40 dex. 
The stars were distributed randomly across the nine 
MOSAIC pointings in Sculptor, in 3 different filters. 
In each image no more than 5% of the total observed 
stars were ever injected as artificial stars at one 
time, so as not to change the crowding properties 
in the image. For each observed pointing 200-400 
images containing artificial stars were created, in 
order to obtain a sufficient total number of artifi- 
cial stars. This resulted in nearly 7000 images con- 
taining a total number of 3.5 million artificial stars 
spread across the full area of the Sculptor dSph. 
Each image was re-reduced using the same tech- 
niques and calibrations as the original image, af- 
ter which the output catalog was matched to the in- 
put catalog (while taking care not to match artificial 
stars to observed real stars), to produce a lookup- 
table containing the recovered artificial stars (and 
the recovered photometric offsets) in 3 filters. 
Figure [3] shows the result of the artificial star tests 
in the I,V-I CMD of Sculptor. The completeness 
level accounts for only those objects that DoPHOT 
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unambiguously recognises as stars. DoPHOT actu- 
ally detects more objects, but the faintest objects 
do not have high enough signal-to-noise to be con- 
fidently distinguished as stars as opposed to unre- 
solved galaxies. Only the unambiguously detected 
stars are considered in the completeness level, since 
those are observed with sufficient accuracy to be 
useful in our analysis. Figure [3] shows the depen- 
dence of the completeness level on both colour and 
magnitude. The dependence on colour in particular 
has to be properly taken into account when generat- 
ing realistic synthetic CMDs. 
The need to carry out artificial star tests over the en- 
tire observed field of view is highlighted in Figure[4] 
Here we show the variation of the 50% complete- 
ness level with distance from the centre of Sculptor, 
in the three different filters. The completeness level 
in the B filter is zero for r c //>0.5 degrees, because 
of the lack of B band observations in the outer parts 
of Sculptor. Figure [4] also shows that the centre of 
Sculptor is less complete at a fixed magnitude limit 
than the outer regions. This is due to the increased 
crowding in the central region, which means it is 
harder to determine the shape of the PSF of stars in 
the centre, causing fewer stars to be unambiguously 
detected at a given brightness level, and placing the 
50% completeness at brighter magnitude levels. 
However, the CMDs at differe nt radii fro m the cen- 
tre (see Figure 5 from lde Boer et al.l201 lb show that 
the photometry in the central part of Sculptor actu- 
ally goes deeper than in the outskirts (due to longer 
exposure times), with a well defined shape of the 
MSTOs. Additionally, using the artificial star test 
results to simulate observational conditions makes 
sure that this changing completeness level is prop- 
erly taken into account when determining the SFH. 
Artificial star tests have only been carried out on the 
long, stacked exposures. For those stars saturated in 
the long exposures, observational effects are mod- 
elled by giving them random offsets within the pho- 
tometric error distribution of stars with similar in- 
strumental magnitude and colour. This is justified 
because the recovery fraction for the brightest stars 
is very close to 100%, and the offsets induced by 
crowding are insignificant and signal-to-noise is not 
an issue. 

The process of including observational errors in 
synthetic CMDs in a statistical manner means that 
they can be directly compared to the observed 
CMDs, to obtain the best matching SFH. 



3.2. Constructing MDF models 

In order to include the spectroscopic MDF at the 
same time as the photometry in Talos we use the 
equivalent of a Hess diagram for [Fe/H]. This is 
done by generating synthetic CMDs of different 
stellar populations and binning the stars in metallic - 
ity. In this way it is possible to create synthetic MDF 
models for stellar populations of different [Fe/H], 
age and [a/Fe]. 

The spectroscopic observations only come from a 
fraction of the RGB, which are in turn a small frac- 
tion of stars present in the CMD. To make sure that 
the synthetic CMDs fully sample the RGB, syn- 
thetic MDF populations are generated with an artifi- 
cially increased initial SFR, which is later corrected 
in the MDF models. To correctly reproduce obser- 
vational limits we construct the synthetic MDFs 
using only stars within a similar magnitude range 
to those of the the observed spectroscopic sample. 
From the number of stars observed both photomet- 
rically and spectroscopically in the same magnitude 
range we then calculate the spectroscopic complete- 
ness fraction on the RGB. This fraction is used to 
scale down the synthetic MDF, and thus match the 
spectroscopic completeness. 

The observed spectroscopic uncertainties on the 
metallicites are simulated by considering each indi- 
vidual synthetic star to have a Gaussian profile with 
a width determined by the average observational un- 
certainty on [Fe/H]. The metallicities of these indi- 
vidual stars are then combined to form a synthetic 
MDF which takes into account observed uncertain- 
ties, and can be directly compared to the observed 
MDF. 

3.3. Determining the SFH 

In Section [XT] we described how Talos can be used 
to create accurate synthetic CMDs, which can be 
compared directly to the observed CMDs. We then 
described in Section [3721 how model MDFs are cre- 
ated in Talos, which add additional information that 
can be used to restrict the SFH. Here we now put ev- 
erything together to determine the SFH of a galaxy. 
The technique used to determine the best SFH min- 
imises the difference between observed and syn- 
thetic CMDs through Hess diagrams to produce the 
cl osest match to the observed data, as described 
in Dotohinl l200l . Photometric Hess diagrams are 
constructed by binning CMDs in magnitude and 
colour space (or three magnitude space for Hess 
diagrams using B,V and I filters) and counting the 
number of stars in each bin. The synthetic Hess dia- 
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grams are scaled to contain a number of stars equal 
to the observed Hess diagrams, to make sure that no 
model is given preference during the SFH determi- 
nation based on the total number of stars it contains. 
In Talos, the synthetic MDFs are used in the fit- 
ting procedure in a similar manner as the photo- 
metric Hess diagrams. The spectroscopic MDF is 
treated as an extra observed dimension, for which 
the difference between model and observations also 
needs to be taken into account. A different weight is 
applied to the photometric and spectroscopic com- 
ponents, in order to enhance the importance of the 
spectroscopic information, which contributes less to 
the overall goodness of fit than the photometry, but 
significantly restricts the possible solutions. 
Since the observed data follows a Poisson distri- 
bution, the goodness of fit is expressed as a differ- 
ence parameter between models an d observations ( a 
Poisson equivalent of^ 2 ), given by Dolphin (2002): 

Xpoisson = 2 Y, m i ~ «i + n l ln ~ 

In which w, is the total number of stars in a synthetic 
Hess bin (m, = £ -SFRyxCMD,-,/ in the photomet- 
ric part, and m, = £ SFRyXMDFjj in the spectro- 
scopic part) and n-, the total number of stars in an 
observed Hess bin. By minimising the difference pa- 
rameter we determine the SFH (given by SFR ; ) that 
best matches both the observed CMD as well as the 
observed spectroscopic MDF. 

3.3.1. SFH uncertainties 

Equally as important as determining the best match- 
ing SFH is determining the uncertainties on this 
choice. The two main sources of statistical uncer- 
tainties on the determination of the SFH are related 
to dat a and parameter samplin g of the SFH solu- 
tion (lAparicio & Hi dalgo 2009I) . 
Data sampling is related to how representative a 
model is of the underlying population. In the con- 
text of CMD analysis, constructing the same model 
twice, with different random samplings of the same 
isochrones will result in two slightly different syn- 
thetic CMDs. In order to quantify how this effect 
may impact the inferred SFH without synthesizing 
an entire new set of synthetic CMDs we generate a 
number of different representations of the observed 
CMD. Synthetic CMDs are generated from the SFH 
obtained in the first run-through of the code and ran- 
domly swapping a certain fraction (by default 20%) 
in the observed CMD. In this way we obtain a new 
"observed" CMD which is representative of a CMD 
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Fig. 5: The input and recovered SFHs of a series of 
short (10 Myr) bursts of star formation at different in- 
put age. The black, solid histogram shows the input SFH, 
given the binning adopted. The red histograms show the 
recovered SFH, along with the fit of a Gaussian distribu- 
tion as the dashed line. The mean and variance (cr) of 
the fitted Gaussian distribution are also listed. 



with a different sampling of the underlying popula- 
tion. By generating a number of these CMDs and 
using them with the original models to find the SFH 
we can quantif y the errors on the SFH caused by 
data sampling (Aparicio & Hidalgo 2009). 
The other major source of error is parameter sam- 
pling. This means that the choice of CMD gridding, 
age bins and metallicity bins will affect the final 
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Fig. 6: The MDF of R GB stars in the Sculptor dSph from Hill et al., in prep (HR sample, blue solid histogram) 
and lKirbv et al] d2010n (red dashed histogram). Panel a) shows the complete samples for both MDFs (HR, going down 
to V«18 and Kirby et al. going down to V«20), while panel b) shows the Kirby MDF truncated at V»18 along with 
the full HR sample. Poissonian errorbars are also shown. 



SFH. In order to quantify the effect of this gridding 
on the recovered SFH, different bin sizes and distri- 
butions are used and then compared. In the parame- 
ter space of age and metallicity three different shifts 
are applied: a shift of half a bin size in age and in 
metallicity and also a shift in both age and metallic- 
ity simultaneously. For the photometric binning two 
different bin sizes are used, with the different shifts 
applied to each. For all these different griddings the 
SFH is determined, after which the results are re- 
turned to a common grid. 

The average of all the different solutions is adopted 
as the final SFH, with errorbars determined as the 
standard deviation of the distribution of solutions. 
This technique has been shown to adequately take 
into account the major uncertainties on the recov- 
ered SFH, and give realistic e rrorbars on the final 
SFH (lAnaricio & Hidalgoll2009h . 

3.4. General properties of Talos 

In order to verify the basic operation of Talos, a 
number of tests have been performed, and the results 
are shown in Appendix[A] We checked our ability to 
correctly recover the age and metallicity of a simple 
synthetic stellar population. We also showed that we 
can accurately recover a synthetic SFH of a contin- 
uous period of star formation from the earliest times 
until now. Furthermore, we tested Talos on observa- 
tions of the globular cluster NGC1904. The metal- 
licity we determined is consistent with that obtained 
from spectroscopic observations, showing that Talos 
can accurately recover the age and metallicity of an 
observed simple stellar population. The main effect 
that makes the output solution different from the in- 



put is the depth of the photometric data, and the set 
of parameter griddings adopted to compute the final 
SFH. These result in limits to the age resolution of 
our final SFH. 

In order to understand the limitations of the final 
SFH determined by Talos, it is important to obtain 
an accurate estimate of the age resolution of the so- 
lution. This is determined by recovering the SFH of 
a set of synthetic pop ulations at different input ages, 
as described by (e.g., Hidalg o et al.ll201 it) . 
Three synthetic populations were generated at the 
distance of the Sculptor dSph, with single short 
bursts of star formation (with a duration of 10 Myr) 
at an age of 8.5, 10.5 and 12.5 Gyr. The metallic- 
ity of the bursts was distributed to match the ob- 
served MDF of the Sculptor dSph, as were the ob- 
servational errors. 

The recovered and input SFH for these bursts are 
shown in Figure [5] The recovered SFH is well fit by 
a Gaussian distribution, which shows that the mean 
of the central peak (fj) is recovered at the correct 
age, with a minor shift of =0. 1-0.2 Gyr. The re- 
covered SFH shows that typically »40% of the to- 
tal input SFR is contained within the central bin, 
at all ages. The star formation is spread out over 
the Gaussian distribution, leading to uncertainties in 
the recovered SFH. However, in the case of constant 
star formation (see Section IA.U the recovered star 
formation rates show a more accurate recovery of 
the input SFH. 

From the Gaussian fits to the three bursts in Figure[5] 
we can derive the variance cr, which determines the 
resolution with which the burst is recovered. The 
three bursts are recovered with an age resolution 
of «1.5 Gyr at an age of 12.5 Gyr, «1.3 Gyr at 
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an age of 10.5 Gyr and Gyr at an age of 8.5 
Gyr, which is consistent with a value of 12% of the 
adopted age. 
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Fig. 7: The relative fractions of different metallicities 
present on the RGB at different photometric depths, with 
respect to an MDF going down to M v =-3.2 (V=19.5) on 
the RGB. The relative fraction indicates by which factor a 
population is incorrectly sampled in the CMD. For com- 
parison the level of the HB is V„ B =20.13 (M v =0.43). 



3.5. Simulating the Sculptor dSph 

To use Talos to determine the SFH of the Sculptor 
dSph the general method described in the preceding 
sections needs to be adapted to the specifics of this 
galaxy. 

3.5.1. General setup 

The observed CMD of Sculptor contains both 
stars that belong to Sculptor as well as Milky 
Way foreground and unresolved background galax- 
ies. To avoid these objects influencing the SFH, 
the observed CMD is "cleane d" using the colour- 
colour diagram presented in Ide Boer et al.l d201 ll) . 
Foreground stars with colours that coincide with 
those of Sculptor are assumed not to significantly 
alter the total number of stars in different evolution- 
ary features, as determined by using the Besan^on 
models to predict the numbe r of Milky Way stars in 
the CMD (Robin et al.ll200l . 
Additionally, th e observed CMDs (shown in 
Figure 5 from Ide Boer et al. I l20Tl show that 
Sculptor contains a significant number of BSS stars, 
which we do not attempt to include in the SFH de- 
termination as a younger popu lation, since they ar e 
likely to be genuine BSS stars (Ma pelii et ai]l2009h . 



To avoid the BSS stars influencing the SFH, we 
generate a mock BSS population which is used as 
a fixed background in the SFH determination. The 
mock BSS stars are constructed by using only the 
main sequence part of populations spanning a range 
in metallicity of -2.5<[Fe/H]<-1.0, with ages be- 
tween 3-4 Gyr. This background Hess diagram is 
subsequently scaled so the number of synthetic stars 
in the BSS region matches the observed number. 
The SFH is determined using the CMD without in- 
cluding the RGB region. The number of stars on 
the RGB is low compared to that in the MSTO 
region, and therefore does not contribute strongly 
to the determination of the SFH. Furthermore, the 
RGB feature has been shown t o substantially var y 
from one set to the other (e.g.. iGallart et al]l2005l) . 
Conversely, there is good agreement between the 
theoretical MSTO phase at old ages (>10 Gyr) in 
different isochrone sets. Instead, information from 
the RGB phase is taken into account through the in- 
clusion of the MDF fitting. 

Determining the SFH without including the spec- 
troscopic MDF shows similar overall trends as ob- 
tained when including the MDF. However, the addi- 
tion of the MDF has a result of putting better con- 
straints on less populated region of the CMD, such 
as the metal-rich stars, which would otherwise be 
fit using a lower metallicity and higher age, due to 
the mismatch between the RGB and MSTO in the 
isochrones. 

A single, well defined distance is assumed for the 
Sculptor dSph in the SFH determination, as ob- 
taine d by a number of re liable distance indica- 
tors (Pietrzvriski et al. 2008). The effect of adopting 
a different distance is shifting the main peak of star 
formation to younger ages for a larger distance and 
older ages for a smaller distance. However, the in- 
ternal distribution of the presented SFHs does not 
change significantly when adopting a different dis- 
tance, given the small observed uncertainties on the 
distance. Furthermore, the good agreement between 
different distance determination methods gives con- 
fidence to the single distance adopted in this work. 

3.5.2. Taking into account the spectroscopic 
information 

The spectroscopic information which is used to pro- 
vide extra constraints on the SFH comes from two 
different sources. For r e ;;<0.2 degrees the metal- 
licity information comes from HR spectroscopy 
on the upper RGB, while in the outer parts of 
Sculptor (0.2<r c //<1.7 degrees) it comes only from 
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Fig. 8: The final (a) SFH and (b) Chemical Evolution History of Sculptor, out to a radius of r eH =l degree, obtained by 
fitting all available photometric and spectroscopic information simultaneously. 
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Fig. 9: The (a) SFH and (b) CEH of the centre of Sculptor, determined using three different isochrone sets. The ob- 
tained solution is shown for the Darthmouth Stellar Evolution Database as the (red) solid line, for the BaSTI/Teramo 
isochrones as the (blue) dotted line and for the Yonsei-Yale isochrones as the (green) dashed line. 



Ca ii triplet spectroscopy which also includes fainter 
stars. Can triplet spectroscopy is also available in 
the central region, which shows that metallicities 
from both samples are place d on the same scale, an d 
consistent with each other (Battag iia et al.l l2008b). 
However, because the spectroscopic sample in 
the central region encompasses only the brightest 
aTmag of the RGB, an intrinsic bias may be present 
in the observed central MDF due to the limited sam- 
pling of the stellar populations on the upper RGB. 
To check for the presence of such a bias in our HR 
spectroscopic sample, a comparison is made with a 
different spectroscopic sample taken in t he same re- 
gion, going down to fainter magnitudes ( Kirb v et al] 
2010). Figure [6^ shows a comparison between the 
normalised MDFs of both samples. The MDF of 
the deeper sample clearly shows a larger fraction of 
metal-poor stars than the brighter sample. Figure^ 
shows that the MDF of both samples look the same, 



when the sample of RGB stars is truncated at the 
same magnitude. Figure|6]shows that an MDF deter- 
mined from only the upper RGB stars in Sculptor re- 
sults in a lower relative fraction of metal-poor stars. 
This effect is due to the luminosity function bias on 
the upper RGB, causing an incomplete sampling of 
the metal-poor components for the brightest RGB 
stars. 

In order to investigate the effects of the bias in 
more detail we consider the sampling of the MDF 
at different photometric depths. Using a synthetic 
CMD of a Sculptor-like galaxy for which the metal- 
licity is known for each star, we construct MDFs 
from RGB stars going down to different photomet- 
ric depths. By comparing to the input parameters 
of the synthetic CMD we find that going down 
to M v =-3.2 (V=19.5) fully samples the overall 
MDF. MDFs of different depths are compared to the 
overall MDF to produce the correct relative frac- 
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tion of stars sampled at each metallicity. Figure [7] 
shows the relative fraction of stars sampled on the 
RGB at different metallicities in CMDs of differ- 
ent photometric depth. The effect of depth on the 
RGB is clearly seen, with shallower MDFs under- 
sampling populations with [Fe/H]<-1.9 dex and 
over-sampling more metal-rich components. A sim- 
ilar effect is also seen when applying the same ap- 
proach to the observed, deep Ca n triplet data, with a 
more severe under-sampling of the populations with 
[Fe/H]<-1.9 dex compared to the synthetic CMD 
results. 

It is important not to neglect the metal-poor com- 
ponent in Sculptor, which is well sampled by the 
MSTO region, but not on the RGB. Figure [7] shows 
that no bias is present in the MDF for [Fe/H] >- 1.9 
dex. Therefore, we adopt this as a metallicity cut- 
off, below which the SFH is only constrained using 
photometry. In the outer parts of Sculptor (r e //>0.2 
degrees) no metallicity cut-off is adopted, since the 
Can triplet spectroscopy extends to faint enough 
magnitudes to avoid this bias. 



3.5.3. SFH parameter space 

To set the limits in age and metallicity used in de- 
termining the SFH of the Sculptor dSph we consider 
all information currently available in the literature. 
We first consider the spectroscopic MDF from Ca n 
triplet and HR spectroscopy. The range in [Fe/H] is 
limited at the metal-poor end by the availability of 
the isochrones, which go down to only [Fe/H] «— 2.5 
dex. This is not a problem, because the spectro- 
scopic MDF shows that the majority of the stars 
in S culptor (92% of the to tal) have [Fe/H] > -2. 5 
dex dStarkenbur g et al.l2010h . The metal-rich is lim- 
ited by the absence of stars with [Fe/H] >- 1.0 dex in 
the observed MDF (see Figure |2). A binsize of 0.2 
dex is assumed for [Fe/H], similar to the average un- 
certainty on [Fe/H] in the spectroscopic sample. 
To constrain [cr/Fe] we use the 89 stars in the 
HR sample in the central 25' diameter region 
of Sculptor. Figure [2] shows that a well-defined 
[a/Fe]-[Fe/H] relation from HR spectroscopy is 
present in the centre of Sculptor. We assume that 
there is no change in this relation with radius and 
directly use it to determine the range in [a/Fe] that 
corresponds to the [Fe/H] used in a particular stel- 
lar population (the region defined by the solid black 
lines in Figure 0. In this way we generate a set of 
populations that cover the full range of [Fe/H] and 
[c/Fe] in the spectroscopic observations. 
To restrict our choice of possible ages we notice 



that the distribution of stars in the observed CMDs 
do not allo w any stars ^5 Gy r old to be present 
in Sculptor (Ide Boer et al.112.0 lib . Assuming a maxi- 
mum age of 14 Gyr, for the age of the Universe, the 
range of ages we consider is thus between 5 and 14 
Gyr old, with a bin size of 1 Gyr. 




24.5 ^ ' — 1 — ' 1 1 5 
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Fig. 10: Synthetic V,B-V CMDs of the MSTO region of 
the Sculptor dSph, colour coded by age as inferred from 
the SFH. 



4. Results 

We have described our method and carried out 
tests to show that it works as expected using real 
and synthetic test data (see Section 13.41 ) Now we 
apply Talos to our photometri c and spectroscopic 
data sets of the Sculpto r Sph dde Boer et all 120111: 
Starkenbu rg et alj|2010l Hill et al., in prep). 
To derive the SFH of Sculptor the available pho- 
tometry (excluding the RGB) is fit simultaneously 
with the spectroscopically determined MDF. 
Within r e ;;<0.43 degrees photometry in three fil- 
ters (B,V and I) is used to determine the SFH, while 
further out only the V and I filters can be used. 
For radii r e //<0.2 degrees the spectroscopic MDF 
obtained from HR spectroscopy is used. This MDF 
is used in the SFH fitting for [Fe/H] >- 1.9 dex, 
to avoid the luminosity function bias described in 
Section 13.51 Below this value only the photometry 
is used to constrain the MDF and SFH. For r e //>0.2 
degrees the MDF obtained from Ca n triplet data is 
used, and no metallicity cut is applied. Furthermore, 
we include a static background Hess diagram to 
account for the presence of a BSS population (as 
discussed in Section 13.5b . The observed HB and 
AGB are ignored, since the isochrone set we use 
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Fig. 11: a) The observed, cleaned I,V-I CMD of the Sculptor dSph. Colours represent the number of stars per bin, on 
a logarithmic scale. To allow comparison with the synthetic CMD, the HB feature is removed from the CMD. b) The 
best synthetic I,V-I CMD of Sculptor, obtained from the SFH determination. The HB and AGB evolutionary features 
are not included in the isochrone set used, and so they are not present in the synthetic CMD. c) The difference between 
the observed and synthetic I,V-I CMD of Sculptor, expressed as a fraction of the observed number of stars. 



does not include these evolutionary features. 
The final Star Formation History and Chemical 
Evolution History (CEH) of the entire Sculptor 
dSph (out to T e u «1 degree) are presented in 
Figure [8] The SFH and CEH display the rate of 
star formation at different age and metallicity 
respectively, in units of solar mass per year or dex 
respectively, over the range of each bin. The total 
mass in stars formed in a bin can be obtained by 
multiplying the star formation rate by the age or 
metallicity range of that bin. 

The overall SFH shows that the Sculptor dSph 
is dominated by an ancient (>10 Gyr old), metal 
poor stellar population. A tail of more metal-rich, 
younger stars is seen at a lower SFR down to an age 
of *6-7 Gyr. The shape of the SFH shows that the 
star formation rate declines with age, suggesting a 
single episode of star formation over an extended 
period of time («7 Gyr). 

The SFH is determined usi ng the Darthmout h 
Stellar Evolution Database (iDotter et alJ [2008). 
To determine the effect of using different 
isochrones sets, the SFH of the central part of 
Sculptor (r e //<0.116 deg) w as obtained using 
three different is o chron e sets ( Potter et al.l 2008 : 
Pietrinferni etafl 12004 12009: lYi etalT I200U 
Kim et all 12002). The obtained SFHs displayed in 
Figure [9] show that the obtained solution is very 
similar in all cases, due to the fact that the SFH 
determination is dominated by the MSTO region, 
for which the different isochrone agree at old ages. 
This gives confidence to adopting the Dartmouth 
isochrone library for the determination of the final 
SFH of Sculptor. 

Using the SFH it is possible to simulate CMDs 
of the Sculptor dSph, which can be used to show 
the distribution of age across the CMD in detail. 



Figure \TU\ shows the synthetic V,B-V CMD of the 
Sculptor MSTOs colour coded with age as inferred 
from the SFH. The separation of different ages is 
clearly visible on the MSTO, which highlights the 
necessity of photometry going down to the oldest 
MSTOs to obtain accurate ages. 




[Fe/H] (dex) 



Fig. 12: The best matching simulated MDF (green) solid 
histogram and observed MDF (red) dashed histogram 
from spectroscopic observations. 



4.1. Reliability of the SFH 

The most basic check of the reliability of the recov- 
ered SFH is to compare the synthetic and observed 
CMDs. Figure [TT] shows the observed I.V-I CMD 
of the Sculptor dSph and the synthetic CMD corre- 
sponding to the best matching SFH. Furthermore, 
the difference between both CMDs is shown, ex- 
pressed as a fraction of the observed star counts in 
each bin. 

The total number of stars in the synthetic CMD 
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Fig. 13: The SFH and CEH (obtained by fitting all available photometry and spectroscopy) for each of the 5 annuli of 
Sculptor, with the corresponding radial extent indicated in the plot. For annulus 5 only the I,V-I CMD is used together 
with the spectroscopy to constrain the SFH. 



is consistent with the observed CMD to within a 
few percent, showing that the total mass in stars is 
well matched to the observations. The CMDs are 
in general a good match, over most of the MSTO 
region. The static background BSS is not repro- 
duced with the same colour distribution as the ob- 
served stars, leading to an under-dense region at 
blue colours (V— 1«0.2) in Figure [TTb . However, 
the transition from the BSS into the MSTO region 
is matched well, despite being modelled in a sim- 
ple way (See Section 13.51 ). Figure [Tit shows that 



the colours of the RGB extend to redder colours 
than observed in Sculptor. The RGB colours of 
the adopted isochrone set are slightly too red for 
the metallicity and age determined from the MSTO 
region and Can triplet spectroscopy, which is a 
known problem f or theoretical isochrones (e.g., 
iGallart et al.l2005h . The HB and AGB are obviously 
not reproduced in the synthetic CMD, as they are 
not included in the Dartmouth isochrone set. 
On the RGB a comparison can be made between the 
observed spectroscopic MDF and the MDF obtained 
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from the SFH determination, as shown in Figure [12] 
for the full extent of the Sculptor dSph. We can see 
that the synthetic MDF is mostly consistent with 
the observed MDF within the errors, as expected, 
given that the MDF is used as an input in the SFH 
determination. A population which is not so well 
fit is the tail of metal-poor stars ([Fe/H]<-2.5 dex). 
From the spectroscopic observations it is known that 
these stars are presen t in Sculptor, although in small 
numbers (»8%, see Starkenb urg et al.l |2010l) . The 
models used to fit the SFH do not contain metal- 
licities low enough to include these stars. Despite 
these uncertainties we feel confident that the dif- 
ferences between the observed and synthetic CMDs 
and MDFs are sufficiently low to conclude that the 
fitted SFH constitutes an accurate representation of 
the full range of ages of stellar populations present 
in Sculptor. 
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Fig. 15: The radial distribution of the total stellar mass per 
square degree in the Sculptor dSph, as determined from 
the SFH modelling. This provides the total mass of stars 
formed over the lifetime of Sculptor, 7.8xl0 6 M Q . 
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Fig. 14: The observed (red) dashed histogram and the syn- 
thetic (green) solid histogram MDFs from the SFH anal- 
ysis in the five annuli. The solid (black) line in the upper 
two panels shows the [Fe/H] above which the HR spec- 
troscopic MDF is unbiased, and can be used in the SFH 
determination. 



4.2. Spatial variations in the SFH 

Due to the extensive spatial coverage of both the 
MSTO photometry and RGB spectroscopy we are 
able to determine the SFH over most of the area 



of the Sculptor dSph. Thus, it is possible to deter- 
mine the radial variation in the SFH in 5 annuli (see 
Figure [131 . each containing a similar number of 
stars observed in the V and I filters. For comparison, 
the Sculptor core radius is 0.1 degrees, just within 
the first annulus. The V and I filters were chosen, 
since they are offer the most complete spatial cover- 
age further from the centre. 

The SFH is determined independently in each of 
these annuli in the same way as for the full galaxy. 
The SFH and CEH of the five selected spatial re- 
gions are shown in Figure [13] using all photo- 
metric filters as well as the available spectroscopy. 
In the outermost annulus the SFH is determined 
only using the I.V-I CMD, together with the spec- 
troscopic MDF. Figure [13] shows that the SFH 
changes significantly with distance from the cen- 
tre of Sculptor. The metallicity an d age gradi- 
ent discussed by Ide Boer et alj d201 lb is quantified 
here. The younger, metal-rich populations are found 
mostly in the central region and drop off towards the 
outer parts. 

Figure [14] shows the observed and synthetic MDFs 
resulting from the SFH determination for all five 
annuli in Sculptor. The MDF shown in the upper 
two panels comes from HR spectroscopy, while the 
MDF in the lower three panels comes from Can 
triplet spectroscopy. Figure [14] shows that the syn- 
thetic MDFs resulting from the SFH determination 
are consistent with the observed MDFs within the 
errors. In the upper two panels the synthetic MDF 
lies above the observed MDF for metallicities lower 
than indicated by the black line and below the ob- 
served MDF for higher metallicities. This is similar 
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to the effects of the luminosity function bias seen in 
Figure |7] 
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Fig. 16: The cumulative radial distributions of stellar pop- 
ulations with different ages. 



Using the SFH at different positions it is possi- 
ble to determine the radial distribution of the to- 
tal mass of stars in Sculptor. This is done by inte- 
grating the SFRs over the duration of star forma- 
tion, at different radii from the centre. Figure Q3] 
shows the radial distribution of the total stellar mass, 
obtained from the SFH modelling. The total mass 
in stars formed is highest in the central parts of 
Sculptor (r e //<0.2 degrees) and drops off steeply to- 
wards the outskirts. The total mass in stars formed 
over all radii in Sculptor for our SFH is 7.8xl0 6 
M Q within an elliptical radius of 1 degree or 1.5 
kpc. To determine the core radius of the total mass 
distribution we use a Sersic profile to fit the distri- 
bution shown in Figure Q3] The resulting core ra- 
dius r c =0. 1 1 +0.04 degrees, is consistent within the 
errors with that derived from the observed density 
profile (Battagli aet al.ll2008al) . 
It is also possible to quantify the radial distribution 
of populations with different ages. Figure [T6l shows 
the cumulative radial distributions of stellar popu- 
lations with different ages, as determined from SFH 
modelling. Figure[l6]is another way of showing that 
the oldest population is the most extended compo- 
nent, while younger populations show an increas- 
in g central concentra tion, as qualitatively described 
in lde Boer et al.l (1201 ll) . 

4.3. Resolving bursty star formation 

One crucial aspect of the interpretation of our SFH 
is: can we distinguish between bursty and con- 



tinuous star formation, given our SFH resolution 
and uncertainties? The intrinsic uncertainties on the 
SFH determination (as described in Section [3.3. Il l 
result in a smoothing of the star formation rates, po- 
tentially hiding the presence of bursts of star forma- 
tion. 

To test our ability to distinguish a series of distinct 
bursts of star formation from continuous star forma- 
tion in Sculptor, we carry out a series of experiments 
using the final SFH of the innermost and outermost 
annuli of Sculptor (anulll and annul5 in Figure [T3l. 
The metallicities and star formation rates from the 
final SFH are adopted, and distributed over one or 
more short bursts of star formation. 
For a single narrow (10 Myr) burst of star forma- 
tion at an age of 13.5 Gyr and -2.2<[Fe/H]<-2.0 
dex metallicity range the recovered and input SFH 
and CEH for annulus 5 are shown in Figure [17] The 
recovered SFH matches well with the final adopted 
SFH of the outermost annulus given in Section [4/2] 
This suggests that the star formation in the outer- 
most part of Sculptor can be approximated with a 
single, short burst of star formation. However, the 
SFH and CEH shown in Figure [17] indicate that the 
spread in age and metallicity of the burst must be 
more extended than assumed here. 
A similar analysis for annulus 1 of Sculptor shows 
that the recovered SFH using a single, short burst 
is clearly inconsistent with the adopted SFH of the 
inner annulus. A comparison of the synthetic CMD 
with the observed CMD of annulus 1 shows a bad 
fit to the data, with a value of x 1 =5. 51 (compared to 
;^ 2 =2.94 for the adopted SFH of annuli), making it 
unlikely to be an accurate representation of the ac- 
tual SFH in the centre of Sculptor. Additionally, it 
seems unphysical to assume that the large spread in 
metallicity observed in the centre of Sculptor (1.5 
dex) can be built up during a single short episode of 
star formation. 

We continue by adding a second burst of star for- 
mation to see if this can improve the match to the 
SFH of the innermost annulus of Sculptor. Figure[T8l 
shows the recovered SFH of a synthetic population 
with two narrow lOMyr bursts of star formation at 
different age, with metallicities and star formation 
rates determined from the adopted SFH of annu- 
lus 1 . In each test an old burst was assumed to take 
place at an age of 13.5 Gyr, with an additional burst 
at a younger age. Figures [T8b,b show that the in- 
put SFH is not recovered as two individual bursts 
if the age separation between the two bursts is less 
than 4 Gyr. Only when the youngest burst is placed 
at an age of 8.5 Gyr, that a separation can be seen 
between the two bursts in the SFH, as shown in 
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Fig. 17: The recovered and input (a) SFH and (b) CEH of a synthetic population with a single short (10 Myr) episode 
of star formation at an age of 13.5 Gyr. The solid black histogram shows the input population given the SFH binning, 
while the red histogram shows the recovered SFH. For comparison, the adopted SFH of annulus 5 of Sculptor is shown 
as the black, dashed histogram. 



Figure [T8b. 

The SFHs shown in Figure[18]are in relatively good 
agreement with the final adopted SFH of annulus 1, 
especially when the youngest burst is placed at 9.5 
Gyr (Figure [T8b). The recovered SFH shows fea- 
tures consistent with the adopted SFH at similar 
strengths. However, a comparison of the synthetic 
CMD to the observed CMD of annulus 1 shows that 
the synthetic populations provide a bad fit to the 
observed CMD, with x 1 values of 4.32, 4.10 and 
4.28 respectively for the populations with the young 
burst at 10.5, 9.5 and 8.5 Gyr. The smooth SFH, as 
given in Section l4~2l provides a better representation 
of the observed CMD, which makes it unlikely that 
strong, short bursts of star formation dominated the 
SFH of the Sculptor dSph. 

A synthetic population with three narrow lOMyr 
bursts of star formation, at 13.5, 10.5 and 8.5 Gyr, 
gives a better agreement with the observed SFH. 
Furthermore, the synthetic CMD provides a bet- 
ter fit to the observed CMD of annulus 1, with a 
X 1 value of 3.96, although still not as good as the 
smooth SFH given in Section l4~2l Taking this ap- 
proach further, adding ever more short bursts of star 
formation will provide a better representation of the 
observed CMD, to the point where we can not dis- 
tinguish between bursty and continuous star forma- 
tion. In the end, this is how our model builds up the 
SFH. The limits of the bursty nature of the Sculptor 
SFH cannot be tied down with CMD analysis at 
such old ages as these with better accuracy than 
shown in Figure [18] However, the form of the CMD 
and the metallicities of individual RGB stars would 
be markedly different if star formation had truly oc- 



curred in a few short bursts, and also other evidence 
such as the broad MDF makes it unlikely. 

4.4. The timescale for chemical evolution 

4.4.1. Determining ages of individual stars 

Using the SFH determined for Sculptor it is possi- 
ble to go back to the spectroscopic sample and de- 
termine the ages of the individual RGB stars. For 
each observed star with spectroscopic abundances 
we find all the stars in a synthetic CMD made with 
the SFH of Sculptor, with the same magnitude (in 
all filters) and metallicity within the observed un- 
certainties. These stars are considered to be repre- 
sentative of the age of the observed star. The mean 
age of the matched sample is adopted as the age of 
the observed star, with the corresponding standard 
deviation of the sample as an errorbar. In the case 
of Sculptor around 100 synthetic stars are typically 
available to compute the mean and standard devi- 
ation for each observed star, ensuring enough stars 
are present to determine a reasonable errorbar. 
To demonstrate the advantage of using the SFH re- 
sults to obtain ages for individual RGB stars a com- 
parison is made to the standard simple isochrone fit- 
ting technique, which uses a fine grid of [Fe/H], age 
and [a/Fe] to obtain the age and uncertainty for an 
observed star along an isochrone. Figure [T9h shows 
the Age-Metallicity Relation (AMR) obtained us- 
ing simple isochrone fitting, and Figure [T9b using 
our full SFH approach. The AMR determined using 
isochrone fitting has bigger age uncertainties and al- 
lows younger ages for each star, which are inconsis- 
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Fig. 18: The input and recovered SFH for a series of syn- 
thetic populations with two short (10 Myr) bursts of star 
formation. The oldest takes place at an age of 13.5 Gyr, 
with the second burst at an age of (a) 10.5 Gyr, (b) 9.5 Gyr 
and (c) 8.5 Gyr. The solid black histogram shows the input 
population, given the SFH binning adopted. The solid red 
histogram shows the recovered SFH along with statistical 
errorbars. For comparison, the adopted SFH of annulus 
1 of Sculptor (Figure [T3t is shown as the black dashed 
histogram. 



tent with the SFH (see Fig [SJ. This is because the 
age for simple isochrone fitting is dependent only 



on the colour of the RGB and the measured [Fe/H], 
which allows a wide range in ages for each observed 
star. Conversely, the more constrained AMR com- 
ing from the SFH is because information from the 
entire CMD (and all modelled evolutionary phases) 
is used to constrain the age of a single observed star. 

4.4.2. Chemical evolution of the Sculptor dSph 

With accurate ages assigned to each RGB star in 
the HR spectroscopic sample it is possible to effec- 
tively measure the evolution of particular elements 
and, for the first time, directly determine timescales 
for chemical evolution. Figure l20h shows [Mg/Fe] 
for RGB stars in Sculptor measured from HR spec- 
troscopy, with the age of the individual stars colour 
coded. Figurel20b shows that [Mg/Fe] follows a uni- 
form trend with age, such that the oldest stars are 
more enhanced than the younger. W e can also de- 
termine the age at whic h the "knee" (Tinslev Hl979l: 
iGilmore & Wvsdl 19911) occurs. This marks the time 
at which SNe la start to contribute to the [Fe/Hl 
content of a ga laxy dMatteucci & Brocatol Il990t 
Matteucci 2003]). In Sculptor, the turnover takes 
place between -1.8<[Fe/H]<-2.0 dex, from which 
we determine that the [Fe/H] from SNe la started to 
be produced 10.9+1 Gyr ago. Therefore, the SNe la 
started contributing noticeably to the chemical evo- 
lution of Sculptor approximately 2+1 Gyr after the 
beginning of star formation. This is the first direct 
measure of this timescale, although it has previously 
been inferred from SNe la timescales an d chemical 
evolution models dMannucci et alJl2006l and refer- 
ences therein). 



5. Conclusions 

We have presented the first detailed SFH of the 
Sculptor dSph, using a combination of deep, wide- 
field multi-colour photometry and spectroscopic 
metallicities and abundances. The method used to 
determine the SFH (described in Section[3]l directly 
combines, for the first time, photometry with spec- 
troscopic abundances. 

The SFH (see Figure [8]) shows feat ures sim- 
ilar t o previous rough SFHs (e.g. Da Costal 
[19841: iHurley-Keller et al] fl999i: iDolphinl 1200*21 
iTolstov et al.l 120031) . but resolves the stellar ages 
with much greater accuracy. Additionally, the SFH 
quantifies the age and metallicity uncertainties and 
provides well motivated errors. The SFH and CEH 
have been determined over a large fraction of the 
Sculptor dSph (^80% of the tidal radius), allowing 
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Fig. 19: The Age-Metallicity Relation of stars on the upper RGB in Sculptor, as determined using (a) simple isochrone 
fitting and (b) from our method incorporating the full SFH and MDF information. 
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Fig. 20: a) [Mg/Fe] measurements for the HR spectroscopic sample of RGB stars in the Sculptor dSph (coloured 
filled circles). The colours represent the age in Gyr, derived from the SFH. Stars in the Milky Way are shown for 
comparison (small grey points), b) [Mg/Fe] is plotted directly against age for the same sample. 



us to quantify the radial age and metallicity gradi- 
ents in Sculptor. 

The SFH shows that star formation took place in 
Sculptor for an extended period of time, from 14 
to 7 Gyr ago (see Figure [HJ with a simple, steadily 
decreasing trend. The spatially resolved SFH over 
the whole galaxy (see Figure [T"3l i shows that a ra- 
dial gradient is present in age and metallicity. The 
MDF (right-hand panels of Figure fT3T > shows that 
the metal-poor, old populations are present at all 
radii while the more metal-rich, younger popula- 
tions are found more toward the c entre, consis- 
tent w ith previous qualitative results dde Boer et al.l 

l^oTTh . 

We explored the temporal resolution of our fi- 
nal SFH for the innermost and outermost parts of 
Sculptor under study. We find that the outermost an- 
nulus is roughly consistent with a single episode of 



star formation. In the inner parts of Sculptor it is 
difficult to find a bursty SFH that matches all the 
available data. Therefore, we find no reason to as- 
sume bursty star formation episodes here and prefer 
the overall picture of a single, extended episode of 
star formation. Additionally, the spatially resolved 
SFH and CEH are more consistent with an age gra- 
dient than with two separate populations. 
The SFH and CEH at different radii from the cen- 
tre are consistent with the scenario where Sculptor 
first experienced a single sizeable burst of star for- 
mation at early times, with more metal-rich popu- 
lations forming ever more concentrated towards the 
centre until about 6-7 Gyr ago, when the star for- 
mation activity ceased. 

Likewise, the chemical evolution of the Sculptor 
dSph seems to be straightforward, according to 
Figure [221 which shows that [Mg/Fe] decreases 



18 



TJ.L. de Boer et al.: Star Formation History of the Sculptor dSph 



steadily with time. The simple decline in the 
[Mg/Fe] distribution suggests that the chemical en- 
richment occurs uniformly over the SFH of the 
Sculptor dSph, with a change in slope when the SNe 
la start to contribute. 

The timescale on which SNe la start to con- 
tribute significantly to the chemical evolution of 
Sculptor (approximately 2+1 Gyr after the begin- 
ning of star formation) is comparable within the 
errors to the timescale expected from the theory 
of SNe la timescales, although inconsistent with 
timescales of prompt SNe la explosions (se e e.g. 
iRaiteri et al.lll99dlMatteucci & Recchill200ll) . The 
small number of stars and the large scatter on the 
knee position in Figure l20h means that the turnover 
age is not exactly defined. A larger sample of metal- 
poor stars is needed to obtain better statistics. 
The formation of a well defined, narrow a- 
element vs. metallicity distribution and a tight Age- 
Metallicity relation (see Figure [19b ) is consistent 
with the overall picture that the Sculptor dSph 
formed stars uninterrupted for an extended period 
of time. The temporal evolution of individual stars 
suggests a slow build up over several gigayears. A 
galaxy with a shorter period of star formation is un- 
likely to show such a well-defined extended AMR. 
If the SFH were more bursty, a greater spread in 
[oyFe] would also be expected, instead of a smooth 
trend with age and [Fe/H] . 

As a result, Sculptor can be considered as a good 
benchmark for isolated star formation over an ex- 
tended period of time at the earliest epochs. 
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Appendix A: Tests of the method 

In order to test the ability of Talos to accurately recover a SFH, a number of tests were made on CMDs for 
which the properties are known, as well as on observed data for a well studied globular cluster. 



A. 1 . Synthetic tests 

First of all, a number of simple artificial stellar populations were generated. By recovering the SFH of these 
populations it is possible to test how accurate our method can recover an input SFH. 



input ■ 
Exact binning 
Full solution f 
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Fig. A. 1 : Upper row: (a) Synthetic SFH and (b) CEH of a single burst of star formation without adding any photometric 
errors. The input SFH is the solid black histogram, while the recovered SFH is the green histogram when using the 
same parameter gridding as the input. The red histogram gives the solution when using a set of different parameter 
griddings. Bottom row: A burst of star formation with realistic errors from artificial star tests. The black histogram in 
(c) the SFH and (d) the ceh is the input SFH and the red histograms indicate the recovered SFH using a set of different 
population gridding schemes. A Gaussian fit to the recovered SFH is shown as a black line, with the mean (p.) and 
variance (cr) also listed. 



We determine the ability of Talos to recover the age, metallicity and SFR in a series of synthetic episodes 
of star formation. The stellar population was generated assuming the distance modulus and reddening of 
Sculptor, and using the artificial star tests to obtain realistic photometric uncertainties. In this way realistic 
colour magnitude diagrams can be obtained, of the same quality as the observed CMD of the Sculptor 
dSph. 

As a first check, we apply the SFH fitting method to a single synthetic episode of star formation. The 
Dartmouth isochrones are used to generate a synthetic CMD with a continuous star formation between 9 
and 10 Gyr and - 1. 9< [Fe/H] <- 1.7 dex. First, the episode was generated without any photometric errors 
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and the SFH determined using the exact parameter gridding of [Fe/H] and age as used to generate the 
population (green histogram). The SFH is also determined using a set of different parameter griddings (red 
histogram, see Figure lATh .b) in order to test the effect of using different griddings to obtain the uncertain- 
ties on the SFH (see Section [3.3.1l l. 

Given exactly the same parameter gridding as the input population, the SFH is recovered at the right age and 
metallicity with the correct strength. A limited amount of "bleeding" is observed, due to the uncertainties 
induced by the quality of the data. The effect of using a set of different griddings is more substantial bleed- 
ing of the star formation rate into neighbouring bins in age and metallicity (as seen in the red histograms 
in Figure lATh .b). This bleeding is a direct consequence of the quality of the photometric data and the use 
of different parameter griddings, and determines the SFH resolution (see Section [3.3.U . Since we do not 
have a priori knowledge of the duration and metallicity of an episode of star formation in the observed data, 
the choice of adopting a single population gridding will lead to biases in the results. Therefore, we always 
consider a range of different parameter griddings in obtaining the final SFH. 



Fig. A.2: (a) SFH and (b) CEH of a continuous episode of star formation including realistic errors. The black line 
indicates the input values and the red histogram the recovered values. 



When including realistic photometric errors determined from artificial star tests for the model of an episode 
of star formation, the same bleeding effect is seen. The parameters of the recovered SFH are determined 
by fitting a Gaussian profile, see Figures lATfc .d. The peaks of the SFH and CEH are recovered well within 
the input values ([Fe/H] mefl „=-1.84 dex, Age mf , fl „=9.64 Gyr), but the star formation is distributed over 
more bins, spreading out the star formation episode in time. Figure [ATI shows that «40% of the total star 
formation is typically recovered within the central peak. Due to the quality of the observed photometric 
data (which only just detects the oldest MSTOs) there remains some degeneracy between the burst strength 
and duration, which could be removed by obtaining deeper CMDs that resolve the MSTO with more accu- 
racy. 

Next we consider a more realistic synthetic population which has experienced constant star forma- 
tion (SFR=Kr 4 M /yr) over the metallicity range -2.5<[Fe/H]<-l.l dex between 6 and 15 Gyr ago. To 
take into account the effect of constraints from spectroscopic observations, the solution was determined tak- 
ing into account a synthetic MDF which samples 50% of the RGB stars. The results are given in Figure [A72l 
It can be seen that the input values are correctly recovered for a synthetic population with constant star for- 
mation, within the SFH uncertainties. 



A.2. Globular cluster NGC 1904 

The final test we carry out is to apply our method to real observations of a Globular cluster, which is, within 
our errors, a simple stellar population. During our observing run Galactic globular cluster NGC 1904 was 
also observed in the B,V and I filters. For NGC1904 there have been several photometric and spectroscopic 
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studies (see Table I A. ft . making it a good test of our method. The reduction of these observations, as well 
as the artificial star tests were done in exactly the same way as for Sculptor. 



Table A.l: Adopted properties of NGC1904 



Property 


Value 


Reference 


(m-M) v 


15.45±0.02 


(Ferraro et al. 1992) 


E(B-V) 


0.01 


(Harris 19961 


[Fe/H] 


-1.579+0.069 dex 


(Carretta et al. 2009) 


[tr/Fe] 


0.31 dex 


(Carretta et al. 2010) 



To obtain the SFH for NGC1904, [a/Fe] was chosen as a fixed value constrained by spectroscopic ob- 
servations (see Table lATt . while the age and metallicity were left as free parameters. The best solution, 
using only the available photometric data is given in Figure IA.3I The SFH and CEH show bleeding into 
neighbouring bins, due to measurement errors and the use of different parameter griddings used. This is 
similar to the bleeding effect seen in Figure |AT1 The SFH clearly indicates a very old population, in good 
agreement with typical globular ages. A Gaussian fit to the metallicity distribution gives a mean value (fj.) 
at [Fe/H]=-1.58 dex, with a variance of <x=0.17 dex. This is in good agreement with the spectroscopic 
[Fe/H] given in Table lA.fl This shows that Talos is able to recover the age and metallicity of a real data set 
with all the errors and uncertainties that it implies. 

These experiments show the capability of Talos to recover the age and metallicity of a stellar population, 
as well as the limit of our ability to unambiguously distinguish a burst of star formation from a more 
continuous SFR over a longer time. 



H = -1.58 
= 0.17 



-1.9 -1.7 -1.5 
[Fe/H] (dex) 



Fig. A.3: (a) The SFH and (b) CEH obtained for NGC1904 from photometric information in the B,V and I filters. A 
Gaussian fit to the MDF is also shown as a black, dotted line, with parameters given in the plot. 
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